Absence of vortex condensation in a two dimensional fermionic XY model 
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Motivated by a puzzle in the study of two dimensional lattice Quantum Electrodynamics with 
staggered fermions, we construct a two dimensional fermionic model with a global (7(1) symmetry. 
Our model can be mapped into a model of closed packed dimers and plaquettes. Although the 
model has the same symmetries as the XY model, we show numerically that the model lacks the 
well known Kosterlitz-Thouless phase transition. The model is always in the gapless phase showing 
the absence of a phase with vortex condensation. In other words the low energy physics is described 
by a non-compact (7(1) field theory. We show that by introducing an even number of layers one can 
introduce vortex condensation within the model and thus also induce a KT transition. 

PACS numbers; 



I. MOTIVATION 

Two dimensional lattice Quantum Electrodynamics 
continues to be of interest today as a test bed for ideas 
and algorithms for lattice QCD P, i, Si,!]. In 
this work we focus on the formulation with staggered 
fermions, and refer to it as LQED2 Q. In the continuum 
limit, the theory is expected to describe the two-flavor 
Schwinger model With massless fermions the two- 
flavor Schwinger model contains an SU{2) x SU{2) chiral 
symmetry. Away from the continuum limit, finite lattice 
spacing effects in LQED2 break the chiral symmetry to 
a C7(l) subgroup. In the mean field approximation this 
symmetry is spontaneously broken. However, since in 
two dimensions strong infrared fluctuations forbid spon- 
taneous symmetry breaking, the mean field result is mod- 
ified 8., J3J : Instead the theory develops critical long range 
(gapless) correlations, which can be detected through the 
chiral condensate susceptibility 
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where ipi and i/jj are the staggered fermion fields at the 
site i on a square lattice. In general one expects x ~ 
in the gapless phase and x ~ ^ when the theory 
develops a mass gap, where A and B are constants. In 
the case of a U{1) symmetric theory, the mass gap can be 
generated only due to vortex condensation. In LQED2, 
at strong couplings one finds that Eq. ([T]) holds with 77 = 
0.5 (lo| . On the other hand in the continuum limit, using 
the result in the two-fiavor Schwinger model, one expects 
77 = 1. Since previous studies have shown no evidence of 
a phase transition from strong to weak gauge couplings 
one might conjecture that the theory is always in the 
gapless phase with 77 varying smoothly from 0.5 to 1.0 as a 
function of the gauge coupling. As far as we know this has 
not yet been shown analytically or observed numerically. 

Although at first glance there does not seem to be any- 
thing strange with the above scenario, a closer examina- 
tion reveals an interesting puzzle. Away from the contin- 
uum limit we expect the low energy physics of LQED2 



to be described by a C7(l) symmetric bosonic field the- 
ory since the fermions are confined. Is this a compact 
or a non-compact U{\) field theory? A compact U{1) 
theory will contain vortices, just like the XY model, and 
these vortices can condense. Thus, it can undergo the 
famous Kosterlitz-Thouless (KT) phase transition, due to 
the presence of vortices, from a gapless phase to a massive 
phase. Further, although the susceptibility x is governed 
by Eq. ([T|) in the gapless phase, r] satisfies the constraint 
< 7? < 0.25. Thus 0.5 < 77 < 1, expected in LQED2, 
appears inconsistent with the physics of a compact (7(1) 
field theory. On the other hand it was shown recently 
that when one studies LQED2 at strong couplings in a 
slab geometry with four two-dimensional layers, one finds 
a KT transition separating a gapless phase and a gapped 
phase. Further in the gapless phase Eq. holds and 
indeed one finds < 7/ < 0.25 as expected [lO]. Thus, we 
conclude that LQED2 behaves like a non-compact (7(1) 
theory on a square lattice, but behaves like a compact 
(7(1) theory on a lattice where four two-dimensional lay- 
ers are coupled to each other. How is this possible? The 
motivation behind our work is to shed some light on this 
question by constructing a simpler model that can be eas- 
ily studied numerically and that clearly exhibits all the 
above features. 

Our paper is organized as follows: In Section |TT] we 
construct a new realization of the lattice XY model with 
fermionic composites. This model is easier to study than 
LQED2 but captures the essential physics. In Section 
mil we show that our simpler model, on a square lat- 
tice, is always in the gapless phase and does not con- 
tain the Kosterlitz-Thouless (KT) phase transition. In 
other words the condensation of vortices is eliminated 
completely suggesting the long distance theory is simi- 
lar to non-compact (7(1) field theory. In Section [TVl we 
show that adding one extra two dimensional layer leads 
to a gapped (massive) phase. This means the extra layer 
introduces vortex condensation. With four or more lay- 
ers the theory contains both the gapless and the gapped 
phases separated by the usual KT transition. In Section 
|V]we present our conclusions and suggest some directions 
for the future. 



II. A FERMIONIC XY MODEL 

The action of the conventional lattice XY (bXY) model 
is given by 

< ij > stands for the sum over nearest neighbor sites and 
Zi = eKp(i(j3i) is a complex variable of unit magnitude at 
site i on a two-dimensional L x L square lattice with 
periodic boundary conditions. The partition function of 
the model is given by 

Z = l[d^]expi-^S). (3) 

Clearly the action and the measure of the partition 
function are invariant under the U{1) symmetry z — > 
exp{i9)z. The bXY model is known to contain vortices 
which leads to the existence of two phases as a function 
of T: a gapless phase where vortices are confined and 
a gapped phase where they are free and condense. The 
phase transition separating them is the well known KT 
transition and occurs at some critical temperature Tc- 
For T <Tc for large L one expects 

x = Y.{^'^-4 + ^r''^-'^L^'' (4) 

where rj depends on T . One of the important predictions 
of the KT transition is that -q 0.25 as (T^ - T) ^ 0. 

Here we construct and study a new type of fermionic- 
XY (fXY) model constructed with Grassmann variables. 
The action of this model is given by 

<ij> <ijkl> 

(5) 

where ipi and V'i are two independent Grassmann vari- 
ables on the site i and < ijkl > stands for all plaquettes 
made up of sites i, j, fc, Z in a cyclic order. The partition 
function is given by 

Z = J [dVc^^]exp(-S'i). (6) 

Note that an overall multiplicative constant 1/T, like in 
Eq. ^ , drops out of all observables up to multiplicative 
constants, if /3 is redefined as T(3. The action and the 
measure in eq. ^ are again invariant under the U{1) 
symmetry ip — > exp(iai9)^p and tp exp(i(Ti0)i/> where 
(7i = 1 on all even sites and —1 on all odd sites. Here 
we ask if the low energy theory of the fXY-model resem- 
bles that of the bXY-model. In particular, as a function 
of f3, are there two phases separated by the KT phase 
transition? 

It is possible to integrate out the Grassmann variables 
and rewrite the fXY-model as a statistical mechanics of 
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FIG. 1: A example of a dimer-plaquette configuration of the 
fermionic XY model. 



bond variables (referred to as dimers) b and plaquette 
variables p. In this representation, the partition function 
is given by 

[b,p] <ijkl> 

where b^ij^ = 0, 1 and p<.ijki> =0,1 are the allowed 
values with the constraint that only one of the four 5's 
or p's associated with a given lattice site must be non- 
zero. Physically, this constraint means that each site 
be connected to only one dimer or one plaquette. An 
example of the dimer-plaquette configuration is shown in 
FigUl 

When j3 — Q the fXY-model defined above is exactly 
solvable and was first studied in [11]. It was found that 
the chiral condensate susceptibility 




satisfies eq. (P) with rj = 0.5. Since 77 does not lie be- 
tween and 0.25, the fXY model is clearly different from 
the cXY model. Below we will argue this difference in 
behavior results from elimination of vortex condensation 
completely. 
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FIG. 2; Chiral susceptibility as a function of L for a variety 
of /3 < 7.5. For larger (3 we find it difficult to reach the ther- 
modynamic limit even though the algorithm may be efficient. 

III. ABSENCE OF VORTEX CONDENSATION 

In order to study the fXY-model we have developed a 
directed path algorithm which is a straight forward ex- 
tension of the ideas presented in [l^. Here we focus on 
the results and postpone the discussion of the algorithm 
to another publication. One of the features of the al- 
gorithm is that it allows us to measure x efficiently. In 
Tab.((T| we compare the results of the algorithm with ex- 
act calculations on small lattices for various values of f3. 



Lattice size 




Exact 


Algorithm 


6x6 


0.0 


3.33640... 


3.3366(9) 


4x4 


0.3 


1.43103... 


1.4307(4) 


4x4 


5.5 


0.34615... 


0.3461(2) 


6x6 


10.0 


0.30663... 


0.3068(5) 


6x6 


100.0 


0.02888... 


0.0288(2) 



TABLE I: Comparison of x obtained using the directed path 
algorithm with exact results. The algorithm appears to re- 
main efficient for relatively large values of P on small lattices. 

In Fig.([T]) we plot x as a function of L for various values 
of /3 ranging from to 7.5. We believe the algorithm 
remains stable for these values of f3. We find that the 
behavior of x as a function of L fits well to the form 
AL^~'^, expected in the gapless phase, for all values of 
p. The values of A and 77 obtained from the fits are 
tabulated in Tab.(in|. As Fig.© and the fits in Tab.(|TT| 
clearly show, there is no sign of a phase transition as a 
function of /? on lattices up to L = 1024. The value of 
r] slowly rises with f}. Based on locality of the theory we 
expect that in the — > 00 limit we should get rj = 2. 
Thus, there is no gapped phase in the fXY model as a 
function of /3: Vortices do not condense and the model 
behaves like a non-compact U{1) field theory. 



f3 


A 


V 


X7D0F 


0.0 


0.242(1) 


0.5 


1.3 


1.0 


0.236(1) 


0.776(1) 


0.5 


2.0 


0.256(3) 


1.000(2) 


0.7 


3.0 


0.294(3) 


1.196(2) 


1.3 


5.0 


0.428(8) 


1.517(4) 


0.3 


7.5 


0.54(1) 


1.765(5) 


1.3 



TABLE II: Values of A and r] obtained by fitting x to the 
form AL^-". 



IV. INTRODUCING VORTICES 

How can we introduce vortex condensation in the fXY 
model? This is indeed possible if we consider > 1 two- 
dimensional layers of square lattices, coupled through a 
local interaction. The action of the A^-layered model is 
given by 

SN = Sl-t ^ -0j;l/>,,iVj,i'V'i,i', (9) 
i,<U'> 

where, in addition to the action given in Eq. ([5]), we 
have added a term that couples neighboring layers I and 
I' represented by < W > at each site i. In the dimer- 
plaquette language the extra term introduces additional 
dimers that connect two neighboring layers at a site. We 
will refer to the fermionic model with N layers as fXY^v. 
Thus the model fXYi = fXY is the model we considered 
above which docs not contain vortex condensation. Here 
we will show that fXY2 and fXY4 contain a gapped phase 
unlike the fXYi model. Here we focus at /3 = and study 
the models as a function of t and N. 

Let us first consider the model with two layers {N = 2). 
Figure [3] shows x a function of L for various values of 
t. Table IIIII gives the value of A and 77 obtained from 
a fit of the data to the form AL^^^. The fits begin to 



t 


A 


V 


XVDOF 


0.01 


0.185(1) 


0.252(1) 


2.2 


0.1 


0.261(1) 


0.259(1) 


0.4 


0.25 


0.337(1) 


0.335(1) 


419 



TABLE III: Values of A and i] obtained from fitting the be- 
havior of X as a function of L, in the range 32 < L < 256, to 
the form AL^''' for different values of t hut N = 2. The fit 
fails for t > 0.25. 

fail for t > 0.25, while they are good for t < 0.1. The 
saturation of x for large L at t = 0.5 is consistent with the 
presence of a finite correlation length which implies the 
existence of a gapped phase due to vortex condensation. 
As t decreases the correlation length increases. Is there 
a phase transition at a finite t? Note that the value of 
77 approaches 0.25 as t is lowered to zero. Interestingly, 
as discussed earlier, 77 — 0.25 is a universal result at the 
critical point in a KT phase transition. Thus, it is likely 
that the correlation length diverges only at t — > and 
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FIG. 3: X a function of L for N = 2. The saturation of x 
for large t implies the presence of a finite correlation length. 
The fits are consistent with a diverging correlation length as 
t 0. 



FIG. 4: Chiral susceptibility as a function of L for A'^ = 4. The 
solid lines for t = 0.4, 0.8 and 1.0 are fits to AL^~^ given in the 
text. This power- law fit begins to fails for t = 1.1 and above. 
In particular, a KT transition is expected at f = ~ 1.05 



not at any finite t. The reason for a good fit for t < 0.1 
could be due to the fact that the correlation lengths at 
these values of t are much larger than the lattice sizes 
explored. 

Next we consider TV = 4, which was studied earlier in 
[13 . In Fig. (ED we plot X as a function of L for vari- 
ous values of t. Table IIVI contains the values of A and 
T] obtained from fits to the data as in the N — 2 case. 
The evidence from the fits is consistent with the following 



t 


A 


V 


XVDOF 


0.4 


0.412(2) 


0.146(1) 


1.3 


0.8 


0.402(2) 


0.188(1) 


0.1 


1.0 


0.396(1) 


0.224(1) 


1.1 


1.1 


0.410(1) 


0.259(1) 


2.8 


1.2 


0.487(1) 


0.338(1) 


124 


1.3 


1.004(5) 


0.586(1) 


1500 



TABLE IV: Values of fitting coefficients A and rj obtained 
by fitting x = AL^~'^ for different values of t. The range 
32 < L < 150 was used in the fit. 



scenario: For t < 1.0 the model is in a gapless phase (ab- 
sence of vortex condensation), while for t > lA the model 
is in a gapped phase (existence of vortex condensation). 
The value of r] close to the critical point, tc ~ 1.05(5), 
is again about 0.25, consistent with a KT phase transi- 
tion. Thus, for iV = 4 the model contains a KT phase 
transition. For even iV > 4 we have evidence, from other 
unpublished studies, that the behavior is similar to that 
of iV = 4. 



V. CONCLUSIONS 

In this work we have constructed and studied a 
fermionic XY model. When the model is studied on a 
square lattice it contains no evidence for a gapped phase 
with vortex condensation. For this reason we believe 
it gives a realization of a two dimensional non-compact 
C/(l) field theory. Introducing an even number of layers, 
leads to a gapped phase phase which must be accom- 
panied by vortex condensation. LQED2, away from the 
continuum limit, gives another realization of the long dis- 
tance physics of our fermionic XY model physics. Our 
work sheds light on why the low energy physics of LQED2 
cannot be described by the bXY model even though the 
symmetries are the same. 

There are several questions that may be of further in- 
terest. For example: 

1. Is it possible to characterize a vortex in the dimer- 
plaquette configuration? We conjecture that a vor- 
tex core may be identified by the existence of a site 
such that, on that site all the layers are connected 
by dimcrs to the neighboring layers. An example 
of such a site is shown in Fig. ([5]) for two layers. 

2. Is there some difference about even and odd values 
of N since our definition of the vortex core does 
not work if N is odd. In one layer we have already 
seen lack of vortex condensation. So what happens 
in TV = 3,5,...? 

3. Does the two layered model approach a KT critical 
point as t — > as conjectured here? Can this result 
be understood through analytic means? 



5 




FIG. 5: A dimer confinguration in two layers. The sites with 
dashed dimers connecting the layers is conjectured to be the 
location of the vortex. Note that for all t > the density 
of such sites will be non-zero leading to a condensation of 
vortices [lo| . 



4. What is the effect of /? > in the layered model. 



Hopefully, this will only introduce more disorder in 
the system. This may be worth studying further 
since a new type of order may be established for 
large (i. 

The answers to these questions could shed more light 
on the field theories described by the dimer-plaquette 
models. 
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